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Abstract. We study the thermodynamic properties of the two-component 2 + 1- 
dimensional massive Dirac fermions in an external magnetic field. The broken time- 
reversal symmetry results in the presence of a linear in the magnetic field part of the 
thermodynamic potential, while in the famous problem of Landau diamagnetism the 
leading field dependent term is quadratic in the field. Accordingly, the leading term 
of the explicitly calculated magnetization is anomalous, viz. it is independent of the 
strength of the magnetic field. The Stfeda formula is employed to describe how the 
anomalous magnetization is related to the anomalous Hall effect. 
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1. Introduction and Model 

Since 80s of the 20th century the condensed matter realization of the (2+1) - dimensional 
Dirac fermions DEI attracts attention of researchers. The discovery of graphene in 
2004 PJ was a tremendous breakthrough in the experimental realization of the massless 
Dirac fermions. Now a family of the solid-state Dirac materials which includes both 
the brothers of graphene, viz. made in 2010 silicene, announced in 2014 germanene 
[4], and cousins such as M 0 S 2 , topological insulators, Weyl semimetals, etc. There are 
even artificially designed in 2012 nephews such as molecular graphene [5] and ultracold 
fermionic atoms in optical honeycomb lattice |6]. An exciting feature of the latter 
artificial condensed matter system is that it allows one the experimental realization [7] 
of the Haldane model [2] where the time-reversal symmetry is broken and a quantum 
Hall effect appears as an intrinsic property of the band structure in the absence of an 
external magnetic field. 

As a matter of fact, the low-energy description of these systems is based on a pair 
of the independent effective Dirac Hamiltonians 

74, (k) = fov F (r]k x Ti + k y T 2 ) + A v t 3 , 


( 1 ) 
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where the Pauli matrices r act in the sublattice (pseudospin) space, the two-dimensional 
wave-vector k = (k x , k y ) is counted off the two independent K,, points (with p — ±) 
in the Brillouin zone, and Vf is the Fermi velocity, and A r; is the mass (gap) term. 
Although the spin degree of freedom is neglected for simplicity, it can be easily included 
if necessary. Accordingly, the most general expression for the gap is given by 


( 2 ) 


A r , = A + p A t , 


where the gap A is invariant with respect to the time-reversal symmetry, and ■//A^ is 
not. 

In the simplest case p], A ? = 0, the mass term A originates from a different on-site 
energy on two triangular sublattices that form hexagonal lattice. The fermion doubling 


[ 8 ] guarantees that the ordinary solid-state Dirac materials are always described by pairs 


of the Hamiltonians (P), so that the full Hamiltonian k) = 'H n=+ i( k) © %, ?= _ 1 (k) 
respects time-reversal symmetry (see Ref. [9] for an overview of the discrete symmetries 
for the various mass terms), 


(n © T 0 )?r (k)(n © r 0 ) = u(-k) 


( 3 ) 


that involves operator n swapping 77 = 1 and ■// = —1 valleys. Here To is the unit matrix. 
The inversion symmetry is, however, broken because the sublattices are inequivalent for 
A^O. Obviously, when a separate K, ; point is considered, the time-reversal symmetry 
is always broken. 

A more sophisticated case, A t 7 ^ 0, with broken time-reversal symmetry is realized 
in the Haldane model [2] by including next-nearest-neighbour (nearest-neighbour on the 
same sublattice) hopping and periodic local magnetic flux density with zero total flux 
through the unit cell. The experimental implementation of this model in [7j is highly 
nontrivial. An earlier history of searches of the condensed matter realization of the 
time-reversal anomaly is presented in [ IU] . 

From a theoretical point of view, the specifics of the Haldane model appears to 
be more transparent when one considers the spectrum of the Hamiltonians (JT|) in an 
external magnetic field B = VxA=(0,0,5) applied perpendicular to the plane along 
the positive z axis. Accordingly, the momentum operator pi = —ihdi has to be replaced 
by the covariant momentum pi ^ pi + ~ c ^-i- Here —e < 0 is the electron charge and 
A = (0, Bx, 0) is the vector potential in the Landau gauge. The corresponding Landau 
level energies are 



( 4 ) 


where M n = yj. A^ + 2nvfh\eB\/c. 

In 2 + 1 dimensions there are two inequivalent irreducible 2x2 representations of 
the Dirac algebra labelled by a = 1, 2. In particular, the Hamiltonian ([TJ) corresponds 
to the representation 7 ^ = (r 3 , — A 7 T 2 , irf) with v = 0,1,2. In general, one can relate 
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the sign of the energy of the n = 0 Landau level to the sign of the product of the mass, 
A,,, and the signature, ?/ a , of the Dirac matrices 7 ^. The signature is defined as mi 


Va 


2 tr ba7a7a] 


±1 


( 5 ) 


and is called sometimes “chirality”. Thus, instead of putting the ± sign in the matrix 
7 *, one can instead include this sign in the mass term (see e.g. Refs. BCG])' 

As mentioned before, in the time-reversal symmetric case, Ax = 0, the spectrum 
of the full Hamiltonian TL( k) = 'H v =+ i(k) © 'H v= _i(k) is particle-hole symmetric and 
invariant under B —> —B. The case Ax 7 ^ 0 and A = 0 corresponds to the so called 
Haldane mass that breaks time-reversal symmetry. As the result, the n = 0 “zero¬ 
mode” level breaks both the particle-hole symmetry and invariance under B —7 —B. 
It is crucial that the difference between these cases survives even in the B — 0 limit, 
making possible the quantum Hall effect in the absence of magnetic field. 

The vast majority of the literature (see e.g. Refs. mmm and references therein) 
on the thermodynamic properties of the Dirac fermions including magnetic oscillations 
is devoted to the case of the particle-hole symmetric spectrum. This concerns both the 
papers that consider the reducible 4x4 representation of the Dirac algebra in 2 + 1 
dimensions and present a field theoretical perspective of the problem ra and the works 
devoted to the condensed matter systems HUES] with even number of fermion species. 
Since the time-reversal symmetry is preserved, a finite Hall conductivity is only possible 
in nonzero external magnetic field. Respectively, the first magnetic field dependent term 
of the grand thermodynamic potential D(/i, B) is proportional to B 2 , so that there is 
no net magnetization in the limit B —> 0. 

The purpose of the present work is to study the thermodynamic properties of the 
Dirac fermions with the broken time-reversal symmetry. In other words, we focus on 
the asymmetric spectrum ([4|) associated with the Hamiltonian ([T]) for one fixed value of 
7], viz. we take A^ = 0, and have A, n = A. Having the thermodynamic potential for 
“unphysical” case with a fixed r), it is straightforward to obtain, for example, the result 
for two K± points with the Haldane mass and include the spin degree of freedom. 

The paper is organized as follows. We already began by presenting in Sec. Q] the 
problem and the model Hamiltonian ([T]). In Sec.[2]we consider the grand thermodynamic 
potential. We show that the unboundedness of the Landau level spectrum for the 
negative energies results in the diverging carrier density and discuss the relativistic 
form of the thermodynamic potential that corresponds to the finite carrier imbalance. 
The magnetic field dependent parts of the carrier imbalance and density are obtained 
in Sec. fusing the Euler-Maclaurin formula. The main results of the present paper that 
include thermodynamic potential and magnetization in the low magnetic field regime 
are presented in Sec. |4j Using the Streda formula we overview in Sec. [5]how the obtained 
expressions can be linked to the known results for the Hall conductivity in the anomalous 
and normal cases. In Sec. [ 6 ] the main results of the paper are summarized. 
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The grand thermodynamic potential ca can be written as follows 


OlT,,/) - -T f rlr.D(r) hi (l - r;V) , 


( 6 ) 


where T is the temperature, /i is the chemical potential, and D(e) is the density of 
states (DOS). We have also set the Boltzmann constant k F = 1. The derivative of the 
thermodynamic potential (( 6 ]) with respect to the chemical potential p, 


p(T,fi) = = J deD(e)n F (e ), 


(7) 


determines the density of carriers p in nonrelativistic many-body theory as a function 
of T, F>, and p. Here n F (e) = l/[exp[(e — p)/T) + 1] is the usual Fermi function and V 
is volume (area) of the system which is set to be unit. 

In the absence of scattering from impurities the DOS per unit area for a given value 
of 77 is expressed via the energies of the Landau levels e nv as follows 

\eB\ 


£>o(f) = ^ [<5(e + rjAsgn (eB)) + £[$(e - M„) + 6(e + M n )] 


71=1 


Substituting the DOS D 0 (e) given by Eq. (jHJ) in Eq. (J7J) we obtain 
\eB\ 


P = P 0 + Pn> 1 — 


27 rhc 


OO 


n F (-rjAsign(eB)) + > \n F (M n ) + n F (-M n )\ 


71=1 


( 8 ) 


(9) 


where p 0 an d p n > 1 denote contributions from the lowest, n — 0 , and remaining (±M n , 
n 0) Landau levels, respectively. We observe that Eq. (j9|) diverges because the 
spectrum ([4]) is unbound for the negative energies e nr) . Since the Dirac Hamiltonian (P) 
has to be regarded as an effective model Hamiltonian derived from the tight-binding 
model, the divergence can be removed by using the appropriate band-width cutoff. 

The other way to remove this divergency is to begin with so called relativistic 
thermodynamic potential (see Refs. [131 HZl HE]) 


D rel (T, n) = -T / deD(e) In I 2 cosh 


TT 


( 10 ) 


To clarify the physical meaning of the potential (TTOl) we differentiate it with respect to 
p and obtain the relativistic carrier density 


p re \T, p) = 


dQ re \T , /j) 
dp 


1 

’2 


deD(e) tanh ■ 


P 


2 T 


( 11 ) 


Then for the particle-hole symmetric spectrum, i.e. when the DOS D(e) is an even 
function of the energy e, D(e) = D(—e) one can see that p rel (T, p = 0) = 0. Further, 
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using the identity 1 = 9(e) + 9(—e) with 9(e) being a step function, one can rewrite the 
last equation in the following form 


OO 

P rel = J deD(e) \n F (e)9(e) - [1 - n F (e)]9(-e )]. 

— OO 


( 12 ) 


This shows that p rel corresponds to the relativistic carrier density or the carrier 
imbalance, viz. p rcl = p + — p_, where p + and p_ are the densities of electrons and 
holes, respectively. Accordingly, when the DOS is the even function of energy, the 
relativistic thermodynamic potential (fTUl) can be presented as a sum of the vacuum, 
electron and hole terms [T3] illustrating its physical meaning. The interpretation of p rcl 
for the case of the asymmetric DOS is discussed in detail in na (see also Ref. [T9] for 
a review). 

Substituting the DOS, D 0 (e) given by Eq. (jHJ) in Eq. (fill) , one obtains 


rel = rel , rel _ \ eB \ 

P ~P0 +Pn> 1 - ^ 


tan h /j + ^Asign(eg)+^ 


2 T 


n =1 


fi-M n n + M n 

tanh ——-h tanh ■ 


2 T 


2 T 


■ ( 13 ) 


where pjj 61 and p^>\ denote contributions from the lowest and remaining Landau levels, 
respectively. One can see that, in contrast to Eq. flU]), Eq. (fT3]i converges. Thus the use 
of the relativistic potential (TO makes the calculation free of divergencies, but at the 
end one should prove that the physical results do not change as compared to the results 
that follow directly from the initial potential ()6|) . 

In general, one can relate the grand thermodynamic potential (J6]) and the relativistic 
potential (TO as follows 

1 f°° 

D(T,p) = D rel (T,p)--y de(p — e)D(e). (14) 


When the D(e) is an even function of the energy e, the last term with the integral in 
Eq. (TO reduces to 

/ OO 

deD(e)9(—e). (15) 

-OO 

Then the carrier imbalance p rcl and total carrier density p are related by the expression 
p rel (T, p) = p(T, p) — phf, where phf is the density of particles for a half-filled band 
(in the lower Dirac cone). Note that both p(T, p) and phf are divergent, so that the 
appropriate cutoff has to be introduced. This reflects the already mentioned fact that 
in the continuum held theory there is no lower bound to the Dirac sea of filled electron 
states. 

The thermodynamic potential (fTOl) can be restored from the particle density 
p rel (p,T, B) by integration 

D rel (T, p, B) = - dep rel (T, e, B) + D C (T, B). (16) 

Here D C (T, B) is a constant of integration. Similarly, integrating p(p, T, B) one restores 
the potential ([6]). 
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The magnetization M in the direction perpendicular to the plane is defined in the 
grand canonical ensemble by the derivative with respect to B at fixed chemical potential, 
i.e. 


M(T,p,B) 


dVL(T, p, B) 
dB 


(17) 


Thus one needs only the magnetic field dependent part of the grand thermodynamic 
potentials (TTOl) and OS]) . A special interest for us represents a linear in magnetic field part 
of the thermodynamic potential those presence would imply a nonzero net magnetization 
in the limit of B —> 0. 


3. Calculation of the carrier imbalance and density 


Since we are interested in the weak field regime, there is no need to evaluate the sum 
over Landau levels exactly. The simplest way to extract the magnetic field dependent 
terms is to use the Euler-Maclaurin formula 

(■oo ^ 

F(n ) « / F(x)dx — —F (0) (18) 

Jo 12 

following the seminal paper on Landau diamagnetism [20]. Firstly we calculate the 
carrier imbalance p T ^>i from Eq. (fT3l) . so that 


1 OO 

-,m +£ 


Pnh 


\eB\ 
47 ike 


£ F M (n) 

n= 1 


(19) 


with 

F ve \n) = tanh ^ + tanh - . (20) 

Repeating the arguments of Refs. [20] |T 6 ] one can see that the use of the Euler-Maclaurin 
formula is justified in the weak-field regime, L(B) -C \p\ with L(B) = \JT\eB\Tw 2 F /c 
being the Landau scale. 

One can see that the term with the integral in Eq. (ITSli is independent of the 
magnetic field similarly to the original Landau consideration [20]. Combining together 
the lowest Landau contribution, pjf 1 and term containing (—l/2)F rel (0), we obtain the 
linear in B contribution, p\ el , to p re] : 


pf = 


\eB\ 
An he 


tanh 


p + ? 7 Asign(ei?) 
2 T 


1 u — |A| 

■- tanh- 

2 2 T 


- tanh 
2 2 T 


( 21 ) 


Analyzing all possible cases eB, pA ^ 0, one can simplify the last expression to the form 


rel = eB sgn(r/A) 
87 rhe 


tanh 


h 


IAI 


2 T 


tanh 


p — |A| 
2 T 


( 22 ) 


Note that there is no linear in B term in the analysis of Landau diamagnetism [20l [16]. 
This term would also be absent if one considers the total contribution from both valleys 
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with r) = ±. Similarly to [2DJ [IS], the term with the derivative (-F rel )'(0) produces the 
quadratic in B part, p\f , of p rel : 


cosh 2 cosh 2 M 

The same approach can be applied to the carrier density (TUT) ■ However, as mentioned 
earlier, the “nonrelativistic” analog of the sum (1191) with 

F{n) = 2[n F (M n ) + n F (-M n )] (24) 

is diverging. To make the series convergent, one can introduce the regularization factor, 
exp(— 5M n ) with 5 —>■ +0. Then the use of the Euler-Maclaurin formula (TT8T) is justified. 
Taking into account the identity n F (e) = 1/2[1 + tanh(/i — e)/2T], one can see that the 
field dependent densities coincide, viz. p\ = pp 1 and pu = pff. 

One can also come to the same conclusion that the magnetic field dependent parts 
of the carrier densities p and p rel coincide by using Eq. (TT4l) . Indeed, taking the derivative 
over p, one can see that 

1 f°° 

p rel (T,/i, B) = p(T, p, B) --J deD(e). (25) 

Here, in contrast to Eq. (fl5|) we do not assume the evenness of the DOS D(e). The 
integral in the RHS of Eq. (125|) corresponds to the total number of states in our bands 
which should be independent of the magnetic field [2IJ- 


P\t = 


(eB) 2 v 2 


96nc 2 \ A\T 


4. Magnetic field dependent part of the grand thermodynamic potential 
and magnetization 


Using Eq. (11611 we restore the thermodynamic potentials, viz. 

O rel (T, p, B) = O rel (T, p,B = 0) + Df(T, p, B) + Q)f(T, p, B) + Q C (T, B), (26) 

and 


fl(T, p, B) = O rel (T, p, B) + n c (T, B). 


(27) 


Here Op 1 and Op 1 are, respectively, the linear and quadratic in B parts of the 
thermodynamic potential obtained by integration of Eq. (1221) for pp 1 and Eq. (1231) for 
pp 1 from — oo to p: 


n?\T,p,B) = 


eHsgn(?/A)T 

47T^,C 


In cosh — — hi cosh ^ 


2 T 


2 T 


(28) 


and 


a\t(T,^B) 


(.cBfvl 
487rc 2 | A | 


tanh 


If + 1^1 

2 T 


— tanh 


h- |A| ~ 

2 T 


(29) 


In Eqs. (06|) and (1271) the arbitrary functions Q C (T,B) and fl c (T,B) may depend on T 
and/or B , but independent of p. 
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Then the magnetization M defined by Eq. (TT7T) in the low-field limit reads 


( 30 ) 


where 


Mi = 


dnf 

dB 


e sgn(r]A)T 
A'kHc 


In cosh 


pi + |A| 
2T 


In cosh 


P 


IAI 


2 T 


(31) 


is the anomalous magnetization, 


dfijf e eBhvp 

dB 247 xhc cl AI 


, u + | A | u — | A | 

tanh-tanh ■ 


2 T 


2 T 


(32) 


is the linear in field part, and the function M C (T, B) is independent of p. Our final result 
(l30D with Mi and Mu given by Eqs. (HTTP and ([32D does not depend on a choice of the 
starting thermodynamic potential that may be taken either (HUD or ([U]). The condition 
that at half-filling the anomalous magnetization is absent, viz. M(T, p = 0, B = 0) = 0 
allows to fix the function M C (T, B ) = 0. 

Persisting in the H ^ 0 limit anomalous magnetization Mi was derived in [22] 
using the low-field expansion for the Green’s function in an external magnetic field. 
This powerful and more complicated method allows one to consider interacting systems, 
while the presented here method uses the explicit form of the noninteracting spectrum 
(|4l). In the T = 0 limit Eq. (HTTP acquires the form 

M i = fa A sgn (p)6(\p\ - |A|) +p sgn (pA)0(| A| - fa|)]. (33) 

For A t = 0 and A v = A, the anomalous Mi contribution from K± points has the 
opposite sign. Accordingly, when the integral contribution of both valleys is considered, 
the leading contribution to the magnetization is linear in B term given by Eq. (132D . In 
the T = 0 limit it reduces to 

Mu = “ / ' 2) ’ < 34 > 


This result is in agreement with the papers |23l E] where the gapped graphene was 
considered, if one takes into account the factor of 4 from the spin-valley degeneracy. 
In the limit A —)■ 0 the corresponding magnetic susceptibility was firstly considered 
by McClure [2Tj in the framework of the studies of diamagnetism of graphite. In 
the opposite limit, L(B) |A|, the vacuum term results in the magnetization 
M(/i = T = A = 0) ~ —\/~B (see [25] [26, |27J E] and references therein). Finally 
we note that the impact of disorder on the magnetization and susceptibility of graphene 
with A = 0 was studied in [28j. 


5. Hall conductivity 

It is shown in [22] that the anomalous magnetization (l3Tf is crucial in obtaining the 
off-diagonal thermal transport coefficient. Below we illustrate how this magnetization 
can be used to obtain the anomalous Hall conductivity and show how to obtain the Hall 
conductivity quantization in graphene from Eq. (HUD . 
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When the chemical potential /j falls in a gap of the energy spectrum, the Hall 
conductivity can be found from the Streda formula |29j : 


^ = - eC (sB)„ = ' eC (w) B ’ <35) 

where the second equality follows from the Maxwell relation. Then either substituting 
the carrier density (122]) in the first equality of Eq. (]35|) or substituting the anomalous 
magnetization (13TT) in the second equality of Eq. (I35|) . one obtains 


a xy (B = 0) 


e 2 sign (rj A) 
8nh 


tanh 


T + I A| 

TT 


— tanh 


If - |A| 

2 T 


(36) 


We stress that this value of the anomalous Hall conductivity does not depend on a choice 
of the starting thermodynamic potential. Furthermore, because Eq. (136]) corresponds 
to the second derivative of the thermodynamic potential with respect to both B and p, 
the presence of an arbitrary function of T and B in H and in the magnetization (1301) 
does not affect the final result. 

In the limit T = 0 the last equation reduces to 

<Jxy = -^sign (t7A)0(|A| - |/i|). (37) 

It is interesting to note that this result corresponds to the long-wavelength limit of the 
static Hall conductivity, viz. lim g _>. 0 lim w _^ 0 a xy (u>, q ) (see [3D] and the discussion in |[5Tj ). 
The usual Hall conductivity, relevant for the transport, is given by an opposite limit, 
crf'y = lim^o liiUg^o cr xy (uj, q). This anomalous Hall conductivity was studied in 132.133] 
and in the clean limit for 77 = + reads 


e 2 sgn(A) J 1, H < |A|, 

4V* \|A|/R M>|A|. 


We observe that Eq. (]38l) agrees with Eq. (]37l) when the chemical potential is inside the 
gap, |/i| < | A| in accord with the conditions of the Streda formula [29] validity. 

Finally we note that the presence of disorder might regularize ED! the uncertainty 
of the order of limits mentioned below Eq. (I37|) . For \p\ < |A| the quantized value of 
the Hall conductivity (]38l) is robust with respect to disorder, while for \p\ > |A| the 
presence of disorder modifies the result for the clean case (see Refs. [32, 33] and [33j for 
a recent discussion). 


5.2. Quantum Hall effect in graphene 

So far we have considered the anomalous Hall effect by taking the limit B — 0. To 
complete our consideration is instructive to consider the Hall conductivity at finite held. 
Since the magnetization (]30l) is valid in the low-held regime, we come back to Eq. (fT3l) 
for the carrier density, because it is valid for an arbitrary strength of the magnetic held. 
The same expression can be directly obtained from the Green’s function by doing the 
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summation over Matsubara frequencies HUESJES]. It is important that, in contrast to 
the the nonrelativistic case, the Matsubara summation in the relativistic case is done [[36] 
without an additional convergence factor. When this factor is present, the Matsubara 
sum produces the Fermi function np as we had in Eq. ([9]), rather than tanh that is 
present in Eq. (1T3h . Note also that to compare the corresponding to Eq. (1T31) expression 
from Ref. [36], one should rewrite the lowest Landau level contribution pg el as follows 


Po el = 


eB 
4:71 HO 


9(eB) tanh 


p + //A 
2 T 


9(—eB) tanh 


ft — //A 
2 T 


(39) 


where we used that \eB\ = eB sgn(eB) = eB[9(eB) — 9(—eB)]. Using that in the T —* 0 
limit tanh(e — p)/T —> 1 — 2 9{—e + p), one can show that Eq. ([3TT1) at T = 0 acquires 
the form 

Po‘ = ^8gnfaA)0(|A| - \p\) + £^sgn^(|p| - |A|). (40) 

Similarly one can consider the T = 0 limit for p^>i part given by Eq. (TT9l) and obtain 


Pn> 1 


|eR|sgn(/i) 
27T he 


n= 1 


(41) 


One can see that the first term of Eq. (1401) corresponding to the anomalous Hall 
conductivity (1361) is time-reversal symmetry breaking, since it is odd under B —>■ —B. 
On the other hand, the second term of Eq. (j40l) and the whole Eq. (T41T) are even under 
B —>- —B. For p — 0 in the limit A —^ 0 the first term of Eq. (1401) leaves behind the 
sgn(? 7 A). This property is called a sign anomaly. As noted in [f?j, this anomaly is 
removed by a finite density effects. Indeed, for p 0 only the second term of Eq. (140T) 
survives in the limit A —^ 0. 

Adding together p\ j el and we recover the results of Lykken et al. [35] and 
Schakel [36] 


p rel = 


eB 


■sgn (pA)P(| A| - \p\) + 


|eR|sgn (p) 


1 

2 + 


Q p 2 - A 2 )c 
2hvf\eB\ 


,(42) 


47 The 27 The 

where [x] denotes the integer part of x. As discussed above, graphene corresponds to 
the case with A^ = 0, A r/ = A —>■ 0. Taking into account the spin-valley degeneracy 
and using the Streda formula we reproduce the half-integer quantum Hall effect 


&xy 


2e 2 sgn (eR)sgn ( p ) 

7T h 


1 

2 + 


p 2 c 


2hvf\eB\ 1 ^ ^ 

The argument given below Eq. (1251) allows to conclude that the same result also follows 
from the nonrelativistic carrier density. The Hall conductivity (]43l) can also be written 
as a function of the free carrier imbalance, Pfri e , 


&xy 


2e 2 sgn (el?)sgn ( p ) 

7T h 


(1 + 

'nhcpf ee ~ 

V2 

[ 2\eB\ 


(44) 


where we used that in the absence of magnetic field = p 2 sgn (p)/(7rh 2 vf). 

Concluding it is worth to stress that the quantum Hall effect is a disorder-induced 
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phenomenon. Accordingly the presented consideration only illustrates the expected 
quantization of the Hall conductivity, rather than explains it. 

6. Conclusion 

In the present work we demonstrated how the Landau approach used to describe the 
diamagnetism of electron gas can be applied for the 2 + 1-dimensional massive Dirac 
fermions. We derived an explicit expression (l30]i for magnetization in the low-held 
regime. It contains held independent (anomalous) term (j3Tj) and the linear in held part 
(l32]h While there exist many approaches that allow one to obtain these terms, the 
use of the Euler-Maclaurin formula as in the original work Landau m is the simplest 
option. Using the Streda formula it is illustrated how the obtained expressions for the 
carrier imbalance (]22l) (or Eq. for the magnetization) and ()42l) are related to the 
anomalous quantum Hall effect and quantum Hall effect in graphene. 
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